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Abstract. We update cosmological bound on axion model. The contribution from the anharmonic effect and the newly 
introduced initial overshoot correction are considered. We present an explicit formula for the axion relic density in terms 
of the QCD scale Aqcd, the current quark masses m^'s and the Peccei-Quinn scale Fa, including firstly introduced 1.85 factor 
which is from the initial overshoot. 
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INTRODUCTION 

The standard model suffers several fine tuning problems. 
They are all related to the symmetry of the nature and its 
breaking. The gauge hierarchy problem which is one of 
the main motivation of the upcoming experiment, LHC 
has its origin in the gauge symmetry breaking sector 
Another problem, strong CP problem is deeply related to 
a discrete symmetry, CP. The fact that the almost of the 
suggestions to solve them are more or less related to the 
invisible missing matter in the Universe is interesting. 

The strong CP problem can be stated as that a CP-odd 
operator 6tr{GG} which cannot be forbidden is unmea- 
surably small ~ 10^" [13]. The introduction of the ax- 
ion can solve this problem beautifully [4, 7]. The axionic 
solution introduces an irrelevant operator -^tr{GG}. It 
can be thought of as making 9 a dynamical field. If then, 
the Vafa-Witten theorem tells that CP conserving point 
a + 6 = becomes a minimum of the free energy by the 
strong interaction [1]. 

However, the introduction of the axion coupling solely 
cannot solve the strong CP problem. One needs one more 
ingredient, the relaxation mechanism for axion. The stan- 
dard big bang cosmology provides it naturally, so the ax- 
ionic solution to the strong CP problem is intrinsically 
cosmological one. This point was recognized soon after 
the suggestion of the axion [8]. The axion is one of the 
possible candidate of the cold dark matter 

After the early studies of cosmological implication 
of the axion. Turner [9] made a quantitative, numerical 
estimate of the axion energy density with a anharmonic 
correction and more precise axion mass formula. 

All of their work depends on the assumption of the adi- 
abatic process. However, the assumption can be broken 
in the two phases of the evolution, at the time of start- 
ing to roll down and at the time of QCD phase transition 
[16, 10]. 



We considered the effects of these two phases on the 
axion energy density and reanalyzed anharmonic effect. 
In addition, many parameters of QCD physics and the 
cosmology are changed [17, 12]. We finally represent the 
explicit formula in terms of these parameters. 



RELIC DENSITY OF AXION 

Just after the inflation, a very light scalar field like ax- 
ion tends to be frozen at a certain value [6], (<I>) = (j)Q = 
const., which is determined stochastically. The dynam- 
ics of a scalar field in the expanding universe with the 
Hubble parameter H{t) is governed by the equation of 
motion 4> + 3//(f)<i> + y'(4>) = 0. The second term in 
the left hand side of it prevents roll down of the scalar 
field, while second term forces it to move to the mini- 
mum. Therefore, after the inflation, the scalar fields stays 
at their initial value until the forces V" to roll down 
overcome the Hubble friction 3H<P. 

If the scalar potential is harmonic one, it starts rolling 
down and oscillates when 3H becomes small as the mass 
m. In the hydrodynamic description, such a coherently 
oscillating scalar fields behaves as like the pressureless 
cold(r = 0) matter If the interaction with the known par- 
ticles is sufficiently weak and its mass is light enough, it 
can be a cold dark matter(CDM) candidate. Since the in- 
visible axion has the interactions and the mass originat- 
ing from quantum effect, they can be small enough. So 
it's a natural candidate of this kind CDM. 

The potential of the axion field can be parameterized 
by V{d) = -C(r)cos(0), where 9= a/Fa and C(r) 
ainstGeV"*(r/lGeV)^". If the initial miss alignment an- 
gle 01 is small enough, it becomes the harmonic poten- 
tial with the mass m„(r) = \JC(J)lF}. At the temper- 
ature Ti satisfying ma{T\) ~ m{T\) = IMgl'h^ /Mp, 
the axion field starts rolling down. In our parametriza- 



2/(4+«) 

^ ^ GeV. Under the adiabatic 



tion, 7i = 1 -!f- 

assumption, H,ma/ma <SC ma, the number of the axion is 
conserved. So the axion energy density at present can be 



easily calculated as pa ^ ma{Ty)ma{Ti)y^^^j 9 

TEMPERATURE DEPENDENT AXION 
MASS 

To give a prediction, we should know oSnst and n. Al- 
though, we cannot solve it fully, all we have to know is 
the mass in the two different regime in which we can esti- 
mate it numerically or phenomenologically. Two regimes 
are zero temperature hadronic phase and the temperature 
enough higher than Aqcd [3]. 

At the zero temperature axion mass is given by 
y/Z/JT+ZjmnfTc/Fa, where Z = ot„/ot^. At the high 
enough temperature, axion mass can be obtained by the 
perturbative instanton calculus under the dilute gas ap- 
proximation [11]. The instanton density with the size p 
is 

2N 



n{p,T) = mum,msCt,{^pf^,(^-^^ .-8-V^^(a) 



1 Mtt^' 

p-'' V r 

^^{-\X-^{2N+Nf)-\2A(X)[l+\(N-Nf)\} ^ 

where B, = 1.33876 and Cn = 0.097163 for 3. With 
the 3-loop running coupling constant, the numerical re- 
sults are well fitted at the T ~ IGeV by the form [2] 

r,. — ^0^:^10-12/' mutn^ms \( Aqcd V'^^ 
«inst - J-yo X lU \^3.6.103MeV3 ) \ 380 MeV ) ' 

7.44259-0.554492(3^ 

POSSIBLE CORRECTIONS 

Our rough estimation reUes on two assumptions, 
the harmonic limit(0i <C 1) and the adiabatic 
approximation(OTa/ma,7/ <C nia). We considered possi- 
ble effects of them on the axion energy density. 

There is also a possible correction from the QCD 
phase transition. However, it turns out to be negligible, 
if one does not include fine tuning [2, 10]. 



anharmonic correction 

In the previous section, we use the fact that the axion 
number is conserved. However, the axion field is in the 
coherent state which is not the number eigenstate. There- 
fore the concept of the number conservation is not obvi- 
ous except for the case of the non-interacting harmonic 



Umit. The conservation quantity, total number of axion 
is actually an adiabatic invariant of the system under the 
slowly varying potential. In the anharmonic limit, ~ 1, 
there is no interpretation Uke the number conservation. 

The adiabatic invariant is given by the area in the 
phase space swept by one period, § p dq [15]. In our 
case, its effect is appeared as a correction factor [2, 14] 
/i = 1^ Vcos0'-cos0i d6'. Now the axion 
energy density with the anharmonic correction is pa — 
ma{Ty)ma{T,)(g^\dlh{6,). 



initial overshoot correction 

There are two adiabatic conditions we assumed. One is 
for the expansion rate to be much smaller than the axion 
oscillation, <C OTq . Another is for the change rate of the 
mass to be smaller than the oscillation. However, at the 
time of starting to roll down, the previous one cannot be 
satisfied because the expansion rate is comparable with 
the mass, 3H ~ ma at that time. 

Therefore the adiabatic approximation does not work 
in the first few oscillations, and the numerical analysis 
must be done in this period. To include its effect in our 
calculations, we defined the another temperature T2 and 
the new misalignment angle Q2 at the temperature T2, 
such that the adiabatic condition is well satisfied after 
T2. We explicitly checked that the values after the first 
half oscillation are good enough. We obtained T2/T1 = 
f2{0i,n) and 9202 = f-i{di,n) numerically, and they are 
presented in the fig. 

As one can sees, the correction factors, /i and /a, 
becomes larger at the large 9\ . It is because of that the 
anharmonic terms make the field stays longer time at the 
around of the top of the potential. So this initial overshoot 
correction includes other anharmonic correction. To split 
them, it is useful to factor out overall constant as like 
/2 (01 , n) ~ C2/2 (01 , n) and /s (0i , n) ~ C2/3 (0i , n), where 
the function with a tilde goes to 1 as 9i is going to zero. 

With the corrections /i, /2 and /s, the relic density 
becomes 



0fc|/|(0i,n)/i(0iC3/3(0i,n)) 

C2 ' /2 (01, n) 
1.846 X 01^x^(01, n). 



where the function 



F(0i,n)=/^(0i,«)/i(0iC3/3(0i,n))//2+"/'(0i,n) 



(1) 




which is shown in Fig. 1 . As one can see, there is a 1 .846 
enhancement in the axion rehc density. 
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FIGURE 1. The combined correction factor F{di,n) 



RESULTS AND CONCLUSION 

With all of corrections we considered, the relic density 
of the axion is estimated as 

p„ ^ 1.449 X 10-" f ^) F(0„«)(eV)^ 

(2) 

where we used Z = mu/m^i ~ 0.5, nij^o =135.5 MeV, 
=93 MeV and gH..(present)=3.91, /(0i) is the anhar- 
monic correction given in Fig. 1 and y is the entropy in- 
crease ratio. 

The energy fraction of the axion is given by 



r, nine f '"""V'W.« \ 

Lla ~ 0.379 X 

V3 -6 -103 MeV/ 

'efF{ei)\ 1^0.701 y 

-0.733 



-0.092 



/ aqcd y-'^Y F,, \ 

V380MeVy VlOi^GeVy 



1.184-O.OlO.T 



(3) 



where x = (AQCD/380MeV) — 1. Of course, these ex- 
pressions gives the bound given in Fig. 2. 
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